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Abstract 

The rapidity distribution of thermal photons produced in Pb+Pb 
colUsions at CERN-SPS energies is calculated within scaling and three- 
fluid hydrodynamics. It is shown that these scenarios lead to very 
different rapidity spectra. A measurement of the rapidity dependence 
of photon radiation can give cleaner insight into the reaction dynam- 
ics than pion spectra, especially into the rapidity dependence of the 
temperature. 

One of the goals of heavy-ion physics is the search for the so-called quark- 
gluon plasma (QGP), a novel phase of matter where quarks and gluons are 
deconfined. It is expected that such a state can be produced, e.g., when 
ordinary hadronic matter is strongly heated or compressed 

Real and virtual photons are promising probes of the QGP (they 
may serve, e.g., as a thermometer |^) since they do not suffer from strong 
interactions. Therefore, their mean free path is large enough for them to 
leave the plasma volume without further interactions. 

Transverse momentum distributions of thermal photons produced in Pb + 
Pb collisions at 160 AGeV were already presented in refs. [|, ^, 0. Recently, 
also the rapidity distribution of electromagnetic radiation was investigated 
in more detail M, U, ITD[ . In these works, however, the collision dynamics was 



simplified by assuming a scaling hydrodynamics solution for the longitudinal 



motion, Vz = z/t |Tl|. This limits the usefulness of the results of ||, |T0[ to 
future collider experiments. In this letter we present rapidity distributions 
of thermal photons in P6(160 AGeV) + Pb (at vanishing impact parameter, 
6 = 0), calculated within the three- fluid hydrodynamical model, and show 
that the photons are very useful to constrain the reaction dynamics at these 
energies. 

Let us first give a brief introduction into the three-fluid model. For a more 
detailed presentation, we refer the reader to refs. ||^, |T^, |T^. The original one- 



fluid hydrodynamic model assumes instantaneous local thermodynamic 
equilibrium in the moment when the nuclei collide and thus is not appro- 
priate to describe the initial stage of ultrarelativistic collisions, at least for 
ELab > 10 AGeV. This problem is solved here by considering more than one 
fluid |T^. The three- fluid model divides the particles involved in a reaction 
into three separate fluids: the first two fluids correspond to the projectile 
and target nucleons, respectively, and the particles produced during the re- 
action are collected in the third fluid. Local thermodynamic equilibrium is 
maintained only in each fluid separately but not between the fluids. The 
fluids are able to penetrate and decelerate each other during the collision. 
Interactions between the fluids are due to binary collisions of the particles in 
the respective fluids. This allows for a treatment of non-equilibrium effects 
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in the initial stage of the colhsion. 

The equation of state (EOS) of the target and projectile fluids is that 
of an ideal nucleon gas plus compression energies. A linear ansatz for the 
compression energy with a compressibility of 250 MeV and a binding energy 
of 16 MeV is used |l| . 

The EOS of the third fluid is that of an ideal gas of massive tt-, p-, u-, 
and r^-mesons. At Tc = 160 MeV we assume a first order phase transition 
into a QGP. For the (net baryon-free) QGP we use the bag-model EOS for 
(pointlike, massless, and noninteracting) u and d quarks. The bag constant 
is chosen such that the pressures of both phases coincide at T = Tc- 

For comparison, we also perform calculations within one-dimensional scal- 
ing hydrodynamics |Tl|] where we also use the latter EOS. Here, the com- 
pressional stage of the collision is not treated and thus two free parameters, 
the initial temperature Tj and (proper) time Tj, where the scaling expansion 



starts, have to be fixed. We use the values given in ref. [|T^: Tj = 300 MeV, 
Ti = 0.22 fm. At later times, the temperature in this model is given by 

T{r) = Tp^i^\ , (1) 



where j = 1, 2, 3 labels the different phases of matter (QGP, mixed phase, and 
hadron gas) and Cj denotes the sound velocity in the corresponding phase. 



4 



All quantities, except for t/^^ and Ti^\ are determined by the equation of 
state, cf., e.g., ref. 0. 

The thermal photon production rate from an equilibrated, baryon-free 
QGP is given (to first order in a and as) by 

where E is the photon energy in the local rest frame of the fluid. In the 
following calculations we fix = g'^/4:7[ = 0.4. As shown in ref. 0, the rate 
for a gas consisting of it-, p-, uj-, and //-mesons may also be parametrized 
by eq. (0). Other contributions, e.g. from the Ai meson as well as the 



effect of hadronic formfactors , are neglected. They are of the same order 
of magnitude as higher order corrections to eq. (j^), which have also not been 
taken into account. Thus, eq. is applied for all phases. The contribu- 
tions from the first two fluids are neglected since, for the reactions considered 
here, these fluids are cooler. Also, they undergo a rapid longitudinal expan- 
sion and thus cool faster than the third fluid. This approximation must, of 
course, break down at large rapidities, where the temperature of the third 
fluid drops rapidly, as discussed below. As a check, eq. (^) has also been 
applied to the projectile and target fluids. This certainly overestimates their 
contribution because they contain only baryons. It turns out that, indeed, 
their contribution to thermal radiation is negligible up to photon rapidities 



The thermal photon spectrum is obtained by an integration over space- 
time: 

kx dkx dy J d^k 
Fig. 1 shows our results. In the three-fluid model the temperature is strongly 
rapidity dependent and so is the spectrum of thermal photons. As already 
noted in ref. there is much more high transverse- momentum radiation 
at midrapidity than in scaling hydrodynamics, which is mainly due to the 
different cooling law. In Bjorken's original model, the temperature is a func- 
tion of proper time only and independent of (fluid-) rapidity, and thus the 
thermal photons show no rapidity dependence. This remains true even if 
transverse expansion p, |T^ is implemented into scaling hydrodynamics. 

At this point we should comment on the different freeze-out procedures 
in the two models. In the three-fluid model the freeze-out of the third fluid 
is done instantaneously at some center-of-mass time tf when its (average) 
temperature drops below Tf = 100 MeV. On the other hand, in the Bjorken 
model, the freeze-out takes place at some proper time Tf defined by T{Tf) = 
Tf = 100 MeV. Thus, the space-time volume for which the photon spectrum 
(^ is determined, is different in both cases. However, at least as far as "hard" 
photons are concerned (i.e., kx ^ Tj) this difference is irrelevant: the shape 



6 



Pb+Pb (160 AGeV) 



10" 



10 

o 



10"- 



2 



10"' 



1-D Bjorkcii, T -300 McV 
1-D Bjorken, Ti=250 MeV 
3- fluid model 



10" 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

YCMS 



Figure 1: Rapidity distribution of thermal photons (for /ct = 2 GeV) calcu- 
lated within the three-fluid model (dots) and longitudinal scaling hydrody- 



namics. 



7 



of the space-time volume differs only for the very late stage {t ^ tf), when the 
temperature is too low to give a sizeable contribution to the "hard" photon 
yield. 

In contrast to the three-fluid model, effects of finite baryon-chemical po- 
tential (which appears due to the expected large amount of baryon stopping 
[O, pn|) on the temperature of the mesons are not accounted for in the 



Bjorken model as presented so far. A finite fiB will mainly manifest itself 
in lower initial temperatures whereas the ne^ effect of /i^ 7^ on "hard" 
(i.e., ^ Ti) photon radiation from the QGP phase was shown to be small 



We assume that this remains true also in the hadronic phase. Using 
the standard argument that relates the final pion multiplicity to the initial 
entropy through the assumption of entropy conservation during the whole 



expansion stage [|T^ E2|, and the expression 



for the entropy density in the QGP phase, it turns out that values 3 < 
fis/T < 9 diminish Tj by 4 — 21%. For simplicity let us assume that 



the ratio of pions to baryons and thus /i^/T is independent of rapidity. The 
effect on the thermal photon radiation is shown in fig. 1: the dotted curve 
was calculated using an initial temperature of 250 MeV, all other parameters 



^That is, the change of the photon rate with /is at fixed T. 



being the same as before. If /is/T is taken to be rapidity dependent this may 
introduce some rapidity dependence into the Bjorken modeL However, the 
resulting photon rapidity spectrum would still lie between the dotted and 
the full curve in fig. 1. 

In principle it would be possible to introduce a (fluid-) rapidity dependent 
initial temperature into the Bjorken model by fitting the pion rapidity dis- 
tribution for all rapidities [§, |10| . Then, eq. (|I]) is applied in each rapidity 
slice separately, i.e. Tj Tiirf), T{t) T{T,ri), and no longer globally. We 
will, however, not adopt this procedure here since it introduces an infinite 
number of parameters and contradicts the original philosophy of scaling hy- 
drodynamics as a simple model with only two free parameters (Tj, Xj). Also, 
this procedure leads to a violation of the conservation laws expressed by 
the hydrodynamic equations of motion unless a finite (and rapidity depen- 
dent) baryon-chemical potential is introduced at the same time 0. It can be 
shown that for longitudinal scaling expansion (even including a cylindrically 
symmetric transverse expansion) dp{T, fi)/dri = 0, or, equivalently, 

df) df) 

has to hold (the partial derivatives with respect to r] being performed at con- 
stant r), where s and n denote the entropy and particle number, respectively. 
From simple arguments one may expect the local photon density to be 
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Figure 2: Squared pion rapidity density. 



proportional to the square of the local pion density. However, in the three- 
fluid model this proportionality is not maintained in the finally observed 
pion spectra. As one can see in fig. 2, the squared pion multiplicity obvi- 
ously decreases much slower with rapidity. Also, the transverse momentum 
distribution of pions is softer than that of photons 0. This can be read- 
ily explained with the help of fig. 3 which shows the rapidity distribution 
of the temperature of the third fluid. Note that the highest temperatures, 
which dominate "hard" photon and high-mass dilepton production, prevail 
only for short times where the rapidity distribution of produced particles is 
narrow (in the present model). This is due to the fact that the third fluid is 
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Figure 3: The (mean) temperature of the third fluid as a function of fluid- 
rapidity. 

initially produced at midrapidity and broadens in rapidity-space during the 
expansion. On the other hand, the final pion distribution emerges at freeze- 
out where the temperature rapidity distribution is much broader. Therefore, 
the final pion (rapidity) distribution is also much broader (in rapidity space) 
than that of the photons. A measurement of the thermal photon rapidity 
spectrum would help to answer if this picture is correct (at these energies), 
or if, instead, the temperature distribution is broad in rapidity space from 
the very beginning of the expansion (as in scaling hydrodynamics). 

In conclusion, the rapidity distribution of thermal photons in Pb + Pb 
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reactions at the SPS has been studied within three-fluid and scahng hydro- 
dynamics. It has been demonstrated that thermal photons provide a power- 
ful tool to constrain the reaction dynamics, i.e., the time and fluid-rapidity 
dependence of the temperature of the hot and dense reaction zone. 
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